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bJO' Abstract 

^' 

.^ ■ We study the quantum phase transition of a A^ two-level system ensemble interacting with an 

^P . optical degenerate parametric process, which can be described by the finite size Dicke Hamiltonian 

plus counter-rotating and quadratic field terms. Analytical closed forms of the critical coupling 

Q^' value and their corresponding separable ground states are derived in the weak and strong cou- 

P3 ' pling regimes. The existence of bipartite entanglement between the two-level-system ensemble 

^ I and photon field as well as between ensemble components for moderate coupling is shown through 

numerical analysis. Given a finite size, our results also indicate the co-existence of squeezed fields 

^ ■ and squeezed atomic ensembles. 
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I. INTRODUCTION 

The study of light-matter interaction has been the central topic of quantum optics; it 
has laid the foundation for laser theory, quantum state engineering, fundamental testing of 
quantum mechanics, and implementation of quantum information processing [l|. Among 
the various systems involving the interaction of photons and atoms, the simplest and the 
most important building block to illustrate interesting quantum phenomena involves just 



a 



the one two-level atom (TLA) [2[ . As the number of TLAs increases, collective effects give 
rise to intriguing many-body phenomena; e. g. the existence of a coherent super-radiant 
phase at zero temperature [3|. 

The Hamiltonian describing the coupling of a non-interacting atomic ensemble with a 
single quantized electromagnetic field mode is equivalent, via a Power-Zineau transformation, 
to the Hamiltonian of a free particle under a field related potential [l|]. The ground state 
energy of such a system is bounded from below by the atomic ground state and the vacuum 
field, ergo a super-radiant phase transition of the ground state in a charge only system 
coupled to a single field mode is not possible j^, |5|. Within the standard minimal coupling, 
long wavelength and rotating wave approximation, and discarding quadratic terms, the 
interaction of a photon field with an ensemble of TLAs is described by the Dicke Hamiltonian, 

-f^Dickc = fujJpO^a + huJaJz + h—={aJ^ + aU_), (1) 

V A'' 



where the transition energy for each one of the A^ TLAs and the radiation field frequency 
are Ua and Up, in that order. The atomic ensemble operators are defined in terms of the 
Pauh matrices for the j-th TLA as J^ = (1/2) X]i=i ^^ ^^^ J± — Tl,j=i '^± ■ The coupling 
strength for the photon- atom interaction is denoted by A. The impossibility for a quantum 
phase transition (QPT) in the source model translates into a violation of Thomas-Reiche- 
Kuhn (TRK) sum rule for an atom if the relation between coupling strength and quadratic 



E 



parameter for a phase transition in the Dicke model were to occur [6|]. In order to observe 
the effect of radiation-matter coupling known as super-radiant phase transition in charge 
only systems, more realistic models, such as spin magnetic moment, statistics or infinitely 
many electromagnetic field modes, should be taken into account [7(]. 

Despite the fact that a QPT is forbidden in the physical system that originated Dicke 
model, the phenomenon is interesting by itself. First, the existence of a QPT in the Dicke 
Hamiltonian was reported as a series of instabilities of the ground state for a finite size 



ensemble of two-level systems [8|. Then, the existence of a super-radiant thermodynamic 
phase transition was proved for an infinitely large ensemble interacting with a coherent 
boson field at a given temperature |9|-[ll[. Also, in the classical limit, studies indicated 
that this quantum critical phenomenon is associated to quantum chaos and ensemble-field 
entanglement 12Nl4|. Non-trivial scaling exponents at the critical point have been discussed 



on the large ensemble size regime [15[. Bipartite intra-ensemble entanglement, due to finite 
size effects, was demonstrated for the Dicke Hamiltonian 16|]. This theoretical understanding 
has motivated proposals for the realization of Dicke model in systems that might allow 

a super-radiant phase transition^ e.o. open dynamical systems involving semiconductor 

quantum wells or quantum dots (l7|, |l8|, open dynamical cavity-QED systems with neutral 



atoms [l9| and ions (20|, and superconducting quantum devices |2l|, |22|, just to mention a 
few. Recently, it has been shown that a standing-wave laser driven Bose-Einstein condensate 
(BEC) coupled to a high finesse optical cavity, accounting for the center of mass motion, 
realizes the Dicke model and that the superradiant phase corresponds to a periodical self- 
organized phase of the atoms 23|, |24 1 . 

Currently, there is great interest in pursuing quantum phase transitions of light since 
photons interacting with atoms should be much easier to study and probe than electrons 
in condensed matter systems 25|, |26|]. Schemes to realize composite Dicke models have 
been proposed; e.g. by combining photon hopping between identical cavities in the photon- 
blockade regime, Mott-insulator to superfiuid QPT has been demonstrated in the Dicke- 
Bose-Hubbard model for an arbitrary number of two- level atoms (271. More exotic QPTs 
of light have been predicted in a Heisenberg spi n 1/ 2 Hamiltonian |28[, two species Bose- 



iubbard model 



arrays of coupled cavities 30|,|3l|, and dual-species optical-lattice cavity 
32| . These studies have brought the possibility to analyze critical quantum phenomena in 
conventional condensed matter systems by manipulating the interaction between photons 
and atoms. 

As nonlinear optics plays an important role in quantum optics, especially in the generation 
of quantum noise squeezed states 33|, |34], a natural question one may ask is how to associate 
nonlinear quantum processes with the phenomenon of QPT. Sub-Poissonian photon statistics 
of the field state, and momentum squeezing of the atomic state, have been predicted for the 
Dicke model 1^. It also has been proposed that squeezing of the photon field carries 



signatures of the associated quantum critical phenomena in the size-consistent Dicke model 



35|. 



In this work, we study the quantum critical phenomena of A^ two-level systems embedded 
within a nonlinear optical medium, with A^ finite. An optical degenerate parametric down 
conversion (PDC) process, where the nonlinear medium is pumped by a classical field of 
frequency 2ujf and that field is converted into pairs of identical photons of frequency Uf 
each, is considered. The corresponding nonlinear interaction Hamiltonian is given by 

^PDC = M«' + «^'), (2) 

where the nonlinear parameter k = x^^^/3 is defined by the second-order nonlinearity coef- 
ficient x^^'' £^nd the classical amplitude of the pumping field (3. By plugging the degenerate 
PDC Hamiltonian, Eq.(l2]), into the Dicke model, Eq.(IT]), and restoring the counter- rotating 
terms we obtain the following Hamiltonian, 

H = hwro^a + hUaJz + h—;^(a + a^) Jx + hnia + a^)^, (3) 



where we have rescaled the system energy by {—hn) and redefined the photon frequency as 
Uf = Up — 2k. In the literature, the Hamiltonian in Eq. (j3]) is that of a non-interacting 
TLAs ensemble driven by an electromagnetic field with the standard minimal coupling in 
the long wavelength limit. Of course, such a QPT does not exist for charge-only systems 
interacting with a finite number of electromagnetic radiation modes 4t3] but the proposed 
nonlinear optics PDC process may provide a strong nonlinearity for the field, which might 



m 



36 



The theoretical quantum 



be coupled to a feasible realization of Dicke model |17H24 

phase transition of the proposed model, Eq. Q, has been shown in the thermodynamic 

imit, i. e. both ensemble size and volume are considered infinitely large, A^ and V — )■ oo. 



39l | . To the knowledge of the authors, the finite size effect on the QPT of this Hamiltonian 
remains unanswered. 

The purpose of this work is twofold. First, we show analytically that a pair of unitary 
transformations effect a rotating wave approximation equivalence on the Hamiltonian in the 



weak coupling regime, thus the exact critical coupling strength can be calculated [16|. In 
the strong coupling regime, we show how the critical coupling strength can be calculated 
under just a semi-classical field approximation. In both cases, the critical coupling value 
of the finite-size A^ Dicke Hamiltonian plus counter-rotating and quadratic field terms is 
independent of the atomic number and agrees with results derived from the classical limit, 



A^ and V^ — )■ oo. Second, treated as a degenerate parametric process, we numerically demon- 
strate the existence of squeezed fields and squeezed atomic ensembles as well as bipartite 
entanglement between the two-level system ensemble and photon field and among ensemble 
components themselves. 

II. WEAK COUPLING REGIME, A < oja- 

In order to derive a set of critical coupling values for the quantum phase transition, 
we first consider the weak coupling regime, i. e., A ^ Ua- In analogy to the unitary 
squeezed operator for a degenerate PDC, we define the following unitary transformation 
and associated parameter, 

f = e<^"-^'"), r] = K/[2{cuf + 2K)]. (4) 

Under the restriction 77 <S 1, the Hamiltonian in Eq. (jH]) is reduced into 

H = f~^Hf, 

^ fujfo'a + h—^Jz + h—=(a + a))Jx + ^fi^. (5) 



Thus, the original Hamiltonian, Eq. (J3]), is approximated by the well-known Dicke Hamilto- 
nian plus the counter rotating terms with a modified field frequency a)/, a modified coupling 
constant g, and a constant energy shift hk, defined as 

u:,= h^^^ \uj, ^=f ^^ + M A, k=( ^^ \n. (6) 

As the total excitation number, N = a) a + J^, does not commute with this Hamiltonian, 
Eq.®, in order to further simplify our problem, a second unitary transformation is used, 

U = e-<^'^''^'y , i = g/[^/N{ua + Qs)l (7) 

where the newly defined parameter fulfills ^ ^ 1 due to the weak coupling regime require- 
ment A ^ uja- Neglecting all but linear powers of the parameter ^, the Hamiltonian in Eq. 
([3]) is written 



w 

A 

N 



hwfa)a + h \uJa + Wa(a + a^)^] Jz + ^~/= [(^■J+ + o.'^J^ j , (8) 



with the extra frequency Ua and the modified couphng A given by the expressions 

^^2 2Qfg 



OJr, 



A = ^^^^^ (9) 



N{uJa + W/) (Wa + iOf) 

The weak regime assumption makes it possible to neglect the extra frequency Ua- Thus the 
weak limit Hamiltonian, Eq.(|H]), is further reduced to the well-known finite size Dicke model 
in Eq. ([1]), of which the ground state can be found exactly and undergoes a phase transition 
at the critical value A = ^ujaUJf 



16 



40|. In our case, the critical coupling value in the weak 



coupling regime can be explicitly expressed as 



Uf{uJa + Uf) + 2K{uJa + 2uJ f) UgjUf + 2k) 

For coupling values fulfilling the condition A ^ \w, we can write the ground state and the 
corresponding energy of our system, 

\Gw) = r|o) (g) 1^), ^ ^J=^ (g) \9)„ 

Ea^ = h[k-^y (11) 

where \g)j denotes the ground state for the j-th TLAs. The ground state is a pure separable 
state and independent on the size of the atomic ensemble. Unlike the Dicke model, here the 
field is in a superposition of vacuum |0) and two-photon |2) states due to the degenerate 
parametric process. 

III. STRONG COUPLING REGIME, A > Ua- 

In the weak coupling regime, the ground state is well described by a finite superposition 
of Fock states times the two-level system ground state; while in the strong coupling regime, 
it is possible to consider the field in a coherent state and try to find the corresponding 
ensemble state. By substituting the photon creation and annihilation operators by their 
expectation values, the Hamiltonian in Eq.(|3]) becomes 

2hX 
Hs = hwf\a\^ + hKa\ + hwaJz + ^=aRJx-, (12) 

V A'' 



where the complex coherent state parameter is defined as a = an. + loti- It is possible to 
arrange this semi-classical Hamiltonian in Eq. flT2|) as a nested array of tensor products of 

6 



the form 



Hs = h [oj f\a\'^ Ka^j^ l2iv + 



I2 + I 



2N-2 



Ho 



h + 1 



2N-1 



H, 



(13) 



where the symbol I^ represents the unit matrix of dimension d and the auxihary matrix of 
dimension two is 

2\aR 



H, = hi^a 



■(Tt. 



N 



(14) 



Thus, the ground state energy is found 



Egs=^ 



N 



uj\a\-' + naj, - -^ul + 16\'al/N 



(15) 



In order to calculate the critical coupling value, we optimize this ground state energy for 
the real and imaginary parts of the coherent state parameter, a, and find the following 
self-consistency equations. 



a 



N[16X^-ul{ujf + AK) 



^ 4A2(w^ + 4fi:)2 

aj = 0. 



(16) 



The phase transition in the strong coupling regime occurs at the critical value given by the 
expression, 



^5 = -^JuJa{uJf+4:K). 



(17) 



Although a finite size has been assumed for the atomic ensemble, this critical coupling value 
found in the strong coupling regime, Eq. 017p . is in accord with that derived from the free 
energy by using the thermodynamic limit method for an infinitely large two-level system 



11| for the reason that in both cases the field is assumed to be in a coherent 



ensemble l6l. llOl 
state. 

The mean-field constrain set, Eq. (TT6|) . approximates, in the strong coupling regime, A 3> 
As, the following ground state and ground state energy. 



N 



\Gs) 



E, 



Gs 



l")QS)l^)i' 

16A2 (tof + 4k) 
7 



(18) 



where the auxihary two-level state is defined as 

|^)= / {\g) + l3\e)), 
^ ;^.(^/ + 4^)-4A^ (19) 

Again, as expected, the ground state is a pure separable state; here the difference is that 
each component of the ensemble is in a superposition of the ground, \g), and excited states, 
|e). Furthermore, for a coupling parameter larger than the nonlinear parameter, A ^ k, 
the auxiliary state is the balanced superposition \v) = {\g) — |e))/-\/2 with null population 
difference, (a^) = 0. 

IV. MODERATE COUPLING REGIME 

Besides the weak and strong coupling regimes, where the ground states are both separable 
states, we apply a direct numerical calculation to find the ground state in the moderate 
coupling regime. In the simulation, each and every single two-level system is taken to be 
indistinguishable from each other and the angular momentum eigenstates basis is used. 



J^\N/2,m) =m\N/2,m), 




J±\N/2,m) = ^/— ( — + 1 ) -m{m±l)\N/2,m±l), 

(20) 

where the Dicke state \N/2,m) is the superposition of all possible ensemble states with 
N/2 + m two-level systems in the excited state and the rest, N/2 — m, in the ground state, 
such that m = -N/2, -N/2 + 1, . . . , N/2 - 1, N/2. 

As the eigenstate of a truncated version of the studied Hamiltonian, Eq.Q, can be easily 
verified to be, or not, an eigenstate for the exact full Hamiltonian, the numerical approach 
taken here consists on assessing a maximum number of allowed excitations for the field, n, 
set to deliver at most a maximum error parameter, e = \Ei — {H)j\/\{H)j\, for a wide range 
of the phase space set by the coupling and nonlinear parameters, (A, k) in that order. The 
set {{Ej, \ipj))} are the numerical eigenvalues and eigenstates, respectively, of the truncated 
Hamiltonian sorted in ascending order, Ej < -Ej+i for j = 1, . . . , (A^ + l)(n -|- 1), and the 
notation {■)j = {ipj\ ■ {ipj) is used. In addition, a degeneracy parameter e^ = \En — Ei\/\Ei\ 

8 



is established to discriminate between non-degenerate and degenerate ground states. In 
the latter case, the proper ground state is constructed as the normalized direct sum of the 
degenerate eigenstates. 

In the following numerical analysis, the error and degeneracy parameter are set to the 
values e < 10~^° and e^ < 10~^°. A maximum of two hundred excitations for the field, 
n = 200, is set in accordance. Numerical results for the on-resonance, Uf = Ua, and off- 
resonance, Uf G [0.85, l)u}a, case are performed for an assorted collection of parameters, 
A^ G [2,6], \/y/N G [0, 5]aJa, k G [0,5]co'a. For the sake of brevity, only those results 
pertaining a bipartite and pentapartite ensemble are shown in Fig. [T]for N = 2 atoms and 
in Fig. [2] for A^ = 5 atoms, respectively. 

The mean value of the z-component of the angular momentum, (Jz), which will be called 
population difference from now on, is shown in Fig. [U^a). Simulation results reveal that, 
as derived in the weak coupling regime, A ^ Ua, the population difference is minimal, 
(Jz) = —N/2, i.e., each and every two-level system is in its ground state, and independent 
of the nonlinear parameter k. Also, for a sufficiently large coupling, A ^ Ua, along with 
an adequate nonlinear parameter k <^ (4A^ — Ua(jJf)/4:Ua such that A ^ k, the population 
difference is null, (Jz) = 0, which relates to the ground state derived in the strong coupling 
regime, Eq. (TTUj) . under the aforementioned restrictions. Fig. [U^b) shows that the numerical 
mean photon number for the field, (n), is in agreement with the general behavior found in 
the analytical results; i.e., the field is in the vacuum field state, with a small two-photon 
component depending on the strength of the nonlinear parameter, for the weak coupling 
regime and in a coherent state, with mean photon number |ap, for the strong coupling 
regime. 

In order to demonstrate the existence of entanglement for the studied Hamiltonian in a 



moderate coupling regime, we ca^ 



culate the maximum shared bipartite concurrence following 



the entangled web approach (4]|, Fig. [I](c), and the field-ensemble entanglement probed 



through von Neumann entropy of the reduced two-level ensemble, also known as entropy 



of entanglement 42|, Fig. [^d). Non-zero regions for both the bipartite concurrence and 
the entropy of entanglement are found in between the separable states corresponding to the 
weak and strong coupling regimes, approximately delimited by the black lines in the Fig. 
Wic) and (d). It is possible to see that the maximum shared bipartite concurrence locates 
in the upper diagonal region. Fig. [Ij^c), indicating that the entanglement shared between 



the ensemble components occurs due to an approximately equal balance between the linear 
atom-photon and nonlinear photon-photon interactions. Instead, the entropy, which has 
its maximum value below the diagonal region. Fig. [U^d), shows a maximum entanglement 
between the two-level system ensemble and the photon field due to a larger atom-photon 
interaction strength. 

A shortcoming of the numerical approach shows up at this point. The area of zero entropy 
below the entangled phase is inversely proportional to the value of the degeneracy parameter, 
edi mentioned above. Also, the error parameter, e, increases as the nonlinear parameter k 
goes to zero. These shortcomings appear due to the truncation of the Hilbert space for 
solving the eigenvalue problem. When the counter-rotating and diamagnetic like terms are 
neglected, the system is confined to certain finite subspaces and the numerical approach 



does not present this problems 16|. Numerical results might be improved by allowing a 
larger maximum for the maximum excitation of the field, optimizing the code, or effecting 
a customized analytical progressive diagonalization scheme based on those presented in 



references 



43 



44|. 



Nevertheless, the current approach allows the calculation of the states for the field and 
atomic ensemble up to the desired precision. In the second row of Fig. [H we show the 
photon number probability distributions, P{n) = \{n\tpg)\'^, related to the four markers, 
labeled from A to D, along a constant coupling parameter, A = 3.323, represented by the 
solid line in Fig. [H^c). In the absence of nonlinear coupling, k = 0, a Poissonian photon 
number distribution is discovered in Fig. [U^e). By calculating the field quadratures variances 
(AX^) and {AY"^), with the field quadratures defined as X = a^ + a and Y = i{a^ — a) , the 
field squeezing is probed. In this case, n = 0, the field is in a coherent state, as expected; 
i.e., the field quadratures mean values and their uncertainty relation have all a value of 
one. With a small value of the nonlinear coupling strength, k = 0.3, the statistics for 
the photon number distribution becomes sub- Poissonian, shown in Fig. [U^f). The field is 
in a squeezed coherent state as the uncertainty relation for the field quadratures remains 
minimal but the variance (AX^) increases as (AF^) decreases. By increasing the nonlinear 
coupling, K = 2.4, an oscillating photon number distribution is found in Fig. Wig)- Now, the 
quadrature squeezing seems to be reversed and the variance (AY^) is smaller than (AX^) 
and close to a value of one. Also, as the value for the quadratures uncertainty relation 
is more than one, the field is no longer in a coherent state. For further increasing of the 

10 



nonlinear coupling, n = 4.8, the oscillating photon number distribution remains. Fig. [TJ^h), 
the variance {AY"^) is further squeezed, and the field is not a coherent state but shows a 
tendency to become the superposition of the vacuum and two-photon state. 

Besides the photon number probability distributions, the corresponding z-component 
angular momentum probability distributions, P{m) = \{'m\tpg)\'^, are shown in the third row 
of FigJH in the same order related to the four points AtoD along the solid line in Fig. [T](c). 
We calculate the mean values and variances for the three momentum operators, (Jj) and 
(AJf) for i = x,y,z, as well as the uncertainty relation between the population difference 
and the dipole phase, 4(AJf)(A<l>^) > 1 where (A<l>^) = (AJ^)/(J^)^. Again, by increasing 
the nonlinear coupling strength. Fig. [T] (i-1), the atomic state changes from a coherent 
atomic state in the absence of the nonlinear parameter, to a squeezed coherent atomic state 
for a small nonlinear parameter. For a larger nonlinear coupling strength, the squeezed 
atomic states becomes a state where the minimal Dicke state, |A^, — A^/2), predominates. 
Our simulation results indicate the co-existence of squeezed fields and squeezed atomic 
ensembles in the moderate coupling regime. The field and atomic statistics for the points 
discussed above, approximated to three decimals for the sake of space, are shown in Table [H 

As the number of two level systems increases, e.g., N = 5 in Fig. O the maximum bipar- 
tite entanglement shared between ensemble components seems to be inversely proportional 
to the ensemble size and the region of entanglement decreases. In the second and third rows 
of Fig. |2l similar photon and atomic statistics from Poissonian, sub-Poissonian, to oscil- 
lating photon number distributions for the field and from the coherent to squeezed atomic 
ensembles, respectively, are demonstrated along a constant coupling parameter, A = 3.019. 

V. CONCLUSION 

Two phase transitions for the ground state were found for a finite size Dicke Hamiltonian 
plus counter-rotating and quadratic field terms, corresponding to the weak and strong cou- 
pling regimes. The ground states before and after these transitions are analytically found 
to be pure separable states, thus there exists no entanglement in the system, identified from 
each other by both the state of the field and two- level system ensemble; i. e., the superpo- 
sition of the vacuum and two photon field states times all the components of the ensemble 
in the ground state, for couplings lesser than the weak critical coupling, and a non-vacuum 

11 



coherent field state times all the components of the ensemble in a superposition of ground 
and excited states, for couplings larger than the strong critical coupling. 

In between these extremes, the ground state presents both ensemble-field entanglement 
and bipartite entanglement between the ensemble components. Results on ensemble bipar- 
tite entanglement behave as expected, the degree of maximum shared pairwise entanglement 
decreases as the number of entangled pairs in the two- level ensemble increases; i.e., for a 
sufficiently large ensemble, e. g., the infinitely large ensemble considered in the thermody- 
namic limit, the maximum shared bipartite entanglement will tend to zero and there will be 
no intermediate region between the weak and strong regimes. Thus, the phase space region 
for which the ground state of the system is entangled is directly related to the finite size of 
the system. 
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D-l I 7:rr- 0.1 



] 




Figure 1. (Color online) The phase diagram of our finite-size Dicke Hamiltonian in the parameter 
space of the linear photon-atom coupling strength, A, and the nonlinear photon-photon interaction 
strength, k. (a) The mean value for the atomic z-component, {Jz), (b) the average photon number 
for the field, (n), (c): the bipartite concurrence, and (d) the entropy of entanglement, (S), are 
calculated for the case of A^ = 2. The corresponding minima and maxima values for the color 
legend are shown below. The field photon number and atomic angular momentum probability 
distributions along the solid line in (c) are shown in (e-h) and (i-1), ordered according to the 
markers A-D, respectively. 
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Figure 2. (Color online) Same as Fig.([T]), but for the case of A^ = 5. 
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2 TLS, Fig. m.e-\) 
A = 3.323 


5 TLS, Fig. de-l) 
A = 3.019 




A 

K = 


B 

K = 0.3 


C 

K = 2.4 


D 

K = 4.8 


A 
K = 


B 

K = 0.3 


C 

K = 2.4 


D 

K = 4.8 


(n) 


22.078 


4.590 


0.502 


0.678 


45.528 


9.417 


0.731 


0.707 


(An2) 


22.084 


3.155 


1.011 


2.135 


45.545 


6.416 


1.153 


2.162 


(AX2) 


1.000 


0.676 


1.824 


2.941 


1.000 


0.676 


1.645 


2.609 


(Ay2) 


1.000 


1.481 


0.869 


0.355 


1.000 


1.480 


1.087 


0.426 


(AX2^(AF2) 


1.000 


1.000 


1.586 


1.045 


1.000 


1.000 


1.789 


1.112 


{J.) 


-1.000 


0.000 


0.000 


0.000 


-2.499 


2.495 


0.000 


0.000 


{^JD 


0.000 


-0.052 


0.903 


0.756 


0.001 


0.005 


5.630 


3.576 


(Jy) 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


{^jD 


0.5000 


0.500 


0.392 


0.364 


1.250 


1.250 


1.147 


0.719 


il) 


-0.023 


-0.512 


-0.466 


-0.781 


-0.068 


-0.153 


-0.861 


-1.902 


i^Jl) 


0.500 


0.5000 


0.487 


0.269 


1.250 


1.246 


1.232 


0.837 


4(AJf)(Al>2) 


1.000 


1.151x10"^^ 


7.515 xlO^'' 


5.887x10^1 


1.000 


1.000 


7.180x103*^ 


1.682x10^1 



Table I. The field and atomic statistics for the markers A to D in Fig. [T]^c) for N = 2 and in Fig. 
[2]^c) for N = 5 two-level systems. See the text for more details. 
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